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^ ■ Abstract 

^; 

The PauU-Villars regularization scheme is often used for evaluating parton distributions 
within the framework of the chiral quark soliton model with inclusion of the vacuum po- 
larization effects. Its simplest version with a single subtraction term should however be 
taken with some caution, since it does not fully get rid of divergences contained in scalar 
and psuedoscalar quark densities appearing in the soliton equation of motion. To remedy 
this shortcoming, we propose here its natural extention, i.e. the Pauli-Villars regularization 
scheme with multi-subtraction terms. We also carry out a comparative analysis of the Pauli- 
Villars regularization scheme and more popular proper-time one. It turns out that some 
isovector observables like the isovector magnetic moment of the nucleon is rather sensitive to 
the choice of the regularization scheme. In the process of tracing the origin of this sensitivity, 
a noticeable difference of the two regularization scheme is revealed. 



^ Email 
^Email 
^ Email 



kubota@kcrn.phys.sci.osaka-u.ac.jp 
wakaniatu@niiho.rcnp.osaka-u.ac.jp 
watabe@rcnp.osaka-u.ac.jp 



1 Introduction 

The recent calculations of nucleon parton distributions within the chiral quark soliton model 
(CQSM) exclusively utilizes the so-called Pauli-Villars regularization scheme [1-6]. This is to 
be contrasted with the fact that most of the past calculations of the nucleon static observables 
were carried out by using the proper-time regularization scheme [7-8] . There are some reasons 
for it. The first reason is mainly technical. For obtaining parton distributions, one need to 
evaluate nucleon matrix elements of quark bilinear operators which is nonlocal in two space- 
time coordinates. The problem is that we have no unanimous idea about how to generalize 
the proper-time scheme for the reguralization of such unusual quantities. The second but 
more positive reason for using the Pauli-Villars regularization scheme has been advocated by 
Diakonov et al. [1,2]. They emphasize that this regularization scheme preserves certain general 
properties of parton distributions such as positivity, factorization properties, sum rules etc., 
which are easily violated by other regularization schemes like the proper-time one. 

Recently, there was a controversial debate on the stability of soliton solutions in the CQSM 
regularized with the Pauli-Villars subtraction scheme [10,11]. It seems that the problem has 
been settled by now, since stable soliton solutions seem to exist at any rate if the Pauli-Villars 
regularization is applied to the quark seas only, not to the discrete bound state sometimes 
called the valence quark orbital. Unfortunately, this is not the end of the story. In fact, soliton 
solutions of the CQSM with use of the Pauli-Villars regularization scheme were obtained many 
years ago by Doring et al. [12]. (To be more precise, the model used by them is not the CQSM 
but the Nambu-Jona-Lasinio model. In fact, they were forced to impose an ad hoc nonlinear 
constraint for the scalar and pseudoscalar meson fields at the later stage of manipulation. 
Otherwise, they would not have obtained any convergent solutions [13].) The fact that the 
single-subtraction Pauli-Villars scheme cannot regularize the vacuum quark condensate was 
already noticed in an earlier paper [14] as well as in this paper [12]. To remove this divergence, 
which is necessary for obtaining a finite gap equation, Doring et al. propose to add some 
counter terms, which depend on the meson fields, to the original effective action. It is very 
important to recognize that this procedure is not workable within the CQSM, since their 
counter terms reduce to mere constants under the chiral circle condition which we impose 
from the very beginning. Thus, one must conclude that the simplest Pauli-Villars scheme with 
the single-subtraction term is unable to fully get rid of the divergence of the vacuum quark 
condensate at least in the nonlinear model. One should take this fact seriously, because it 
brings about a trouble also in the physics of soliton sector. To understand it, one has only to 
remember the fact that the scalar quark density appearing in the soliton equation of motion is 



expected to approach a finite and nonzero value characterizing the vacuum quark condensate 
as the distance from the sohton center becomes large [15]. This necessarily means that the 
scalar quark density appearing in the soliton equation of motion cannot also be free from 
divergences. 

The purpose of the present study is then twofold. On the one hand, we want to show that 
the single-subtraction Pauli-Villars scheme is not a fully satisfactory regularization scheme, 
and that at least one more subtraction term is necessary for a consistent regularization of 
the effective theory. This will be made convinced through the formal discussion given in II 
and also the explicit numerical results shown in III. A. On the other hand, we also want to 
know the regularization-scheme dependence of the CQSM through the comparative analysis of 
typical static observables of the nucleon predicted by the two regularization schemes, i.e. the 
Pauli-Villars one and the proper-time one. The discussion on this second issue will be given 
in III.B. We then summarize our conclusion in IV. 



2 Pauli-Villars regularization scheme 

We begin with the effective lagrangian of the chiral quark model with an explicit chiral sym- 
metry breaking term as 

C-CQM = Cq^ C , (1) 

where £o denotes the chiral symmetric part [16] given by 

£o = ip{i ^-MU^"{x))iIj, (2) 

with 

while 

C'=\f^mltT{U{x) + U\x)-2), (4) 

is thought to simulate a small deviation from the chiral symmetric limit. Here the trace in (4) 
is to be taken with respect to flavor indices. (One could have taken an alternative choice that 
introduces explicit chiral-symmetry-breaking effects in the form of quark mass term. We did 
not do so, because it turns out that this form of action cannot be regularized consistently with 
the Pauli-Villars subtraction method.) 



The idea of the Pauh-Villars regularization can most easily be understood by examining 
the form of the effective meson action derived from (1) with the help of the standard derivative 
expansion : 

S,ff[U] = Sf[U] + Sm[U], (5) 

where 

Sf[U] = -iN^Sp\og{i0-MU^') 

= I d^x {ANcM'^hiM) ti{d^Ud^U^) + higher derivative terms}, (6) 

S„AU] = j d'x\flmlii{U{x) + U\x)-2). (7) 

In eq.(6), the coefficient 



hiM) = -I 



of the pion kinetic term diverges logarithmically. In fact, by introducing a ultraviolet cutoff 
momentum a that should eventually be made infinity, one finds that 

/2(M)~^{lna'-lnM2-l}. (9) 

This logarithmic divergence can be removed if one introduces a regularized action as follows : 

S:ff[U]=SY'[U]+SrrAU], (10) 

where 

srp]^Sf[u]-[j^j 5f-[f/]. (11) 

Here Sf ^^ is obtained from Sf[U] with M replaced by the Pauli-Villars regulator mass Mpv- 
Further requiring that the above regularized action reproduces correct normalization for the 
pion kinetic term, one obtains the condition : 



In ^:;^ =fl (12) 



which can be used to fix the regulator mass Mpy. Once the effective action is regularized, the 
static soliton energy should be a finite functional of the soliton profile F{r) under the standard 
hedgehog ansatz U{x) = exp[iT ■ rF{r)]. Since the soliton equation of motion is obtained 
from the stationary condition of the static energy against the variation of F{r), everything 
seems to be going well with the above single-subtraction Pauli-Villars regularization procedure. 
Unfortunately, this is not the case. To understand what the problem is, we first recall the 
fact that the scalar quark density appearing in the soliton equation of motion is expected to 
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approach a finite and nonzero constant characterizing the vacuum quark condensate as the 
distance from the sohton center becomes large [15]. (This is a natural consequence of our 
demand that both of the soliton {B = 1) and vacuum(i? = 0) sectors must be described 
by the same (or single) equation of motion.) On the other hand, it has been known that the 
vacuum quark condensate contains quadratic divergences that cannot be removed by the single- 
subtraction Pauli-Villars scheme [12,14]. This then indicates that the scalar quark density 
appearing in the soliton equation of motion cannot also be free from divergences. 

To get rid of all the troublesome divergences, we propose here to increase the number of 
subtraction terms, thereby starting with the following action : 

S:ff[U]=SY'[U]+Sr^[U], (13) 

where 

N 

S}^^[U]^Sf[U]-Y.c^SfnU], (14) 

1=1 

with A^ being the number of subtraction terms. The logarithmic divergence of the original 
action is removed if the condition 

l-g.(jj) ^0 (15) 

is fulfilled. Similarly, the normalization condition (12) is replaced by 

The single-subtraction Pauli-Villars scheme corresponds to taking A^ = l,Ai = Mpy, and 
Ci = (M/Mpv) ■ This is naturally the simplest case that satisfies both conditions (15) and 
(16). 

To derive soliton equation of motion, we must first write down a regularized expression for 
the static soliton energy. Under the hedgehog ansatz n{x) = fj^rF{r) for the background pion 
fields, it is obtained in the form : 

E::!uc[F{r)] = Ey^[F{r)] + E^[F{r)], (17) 

where the meson part is given by 

E^[F{r)] = -flml j d\ {cos F{r) - 1) , (18) 

while the fermion (quark) part is given as 

Ey^[F{r)] = E^^, + E^-^, (19) 



with 

E,ai = N,Eo (20) 

N 

K;' = N^T.{En-E^^)-Y.c.KY.{Et-E^^^^). (21) 

n<0 i=l n<0 

Here En are the quark single-particle energies, given as the eigenvalues of the static Dirac 
hamiltonian in the background pion fields : 

H\n>=En\n>, (22) 

with 

H = — — + (3M (cos F(r) + i^^r ■ r sin F(r)) , (23) 

while the energy E^^^ denote the energy eigenvalues of the vacuum hamiltonian given by eq.(23) 
with F{r) = or f/ = 1. Eq.(19) means that the quark part of the static energy is given as a 
sum of the contribution of the discrete bound-state level and that of the negative energy Dirac 
continuum. The latter part is regularized by subtracting from the Dirac sea contribution a 
linear combination of the corresponding sum evaluated with the regulator mass Aj instead of 
the dynamical quark mass. {E^^ in these subtraction terms are the eigenenergies of the Dirac 
hamiltonian (23) with M replaced by Aj and with the same background pion field.) 

Now the soliton equation of motion is obtained from the stationary condition of El^ll^^^[F (r)] 
with respect to the variation of the profile function F{r) : 
6Estatic[F{r)] 







which gives 



5F{r) 
Awr^ {-M [S{r) sinF(r) - P{r) cosF(r)] + /^m^sinF(r)} , (24) 



F{r) = arctan I ^^ ^^''\,^2 I • (25) 



S(r) - ii^ 



M 

Here S{r) and P{r) are regularized scalar and pseudoscalar densities given as 

^ A- 
S{r) = S.ai{r) + j:Snir)-J2c,-^j:stir), (26) 

n<0 i=l ■''^ n<0 

N A. 

P{r) = P.ai{r) + Y.Pn{r)-Y.^.j^T.Ptir), (27) 

n<0 j=l ^^^ n<0 



with 



— / d^x < n\x > Y -^^—^n 
An J I ' ^2 

Nc f ,-^ , . n -^ H\x\ -r) 

— / d X < n\x > t'j l5 T ■ r ' 
Air J ' r^ 
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Snir) = -^ I d^x < n\x > -f"^ ^^ l^ ^ <x\n>, (28) 

Pn{r) = -j^ I dx < n\x > iYl5T -r ^ -<x\n>, (29) 



and Svai{r) = S'„=o(^) and Pvaiij') = Pn=o{r), while 5'^'(r) and P^'if^) are the corresponding 
densities evaluated with the regulator mass Aj instead of the dynamical quark mass M. As 
usual, a self-consistent soliton solution is obtained in an iterative way. First by assuming an 
appropriate (though arbitrary) soliton profile F{r), the eigenvalue problem of the Dirac hamil- 
tonian is solved. Using the resultant eigenfunctions and their associated eigenenergies, one 
can calculate the regularized scalar and pseudoscalar quark densities S{r) and P{r). Eq.(25) 
can then be used to obtain a new soliton profile F{r). The whole procedure above is repeated 
with this new profile F{r) until the self-consistency is fulfilled. 

Now we recall an important observation made before. The scalar quark density S{r) at 
the spatial infinity r = oo with respect to the soliton center should coincide with the scalar 
quark density in the vacuum {B = 0) sector, which is nothing but the familiar vacuum quark 
condensate (per unit volume) {'4''4')vac- That is, the following simple relation must hold : 

{m.ac = ^ J Sir = 00) dh = S{r = 00). (30) 

(Later, this relation will be checked numerically.) What we must do now is to find necessary 
conditions for the subtraction constants q and Aj in the multi-subtraction Pauli-Villars scheme 
to make the vacuum quark condensate finite. This can be achieved by examining the expression 
of the vacuum quark condensate obtained consistently with the soliton equation of motion : 

N / A \ 

Mi^^YZ = Mii^^l^^ - E Q (^) AS^)t, (31) 

or equivalently 

m:z = {^^).ac - E c. (^) {^^)t, (32) 

where 

with E]^' = (A;2 -|- M^)"^'^, while {i^i^Y^ are obtained from (V^V')t,ac with the replacement of 
M by Aj. Using the integration formula 

_fk_ 1_ 

(27r)^ Vk^ + M2 
with a being a ultraviolet cutoff momentum, we obtain 

Z7I - I ~ — ; \ivi / ~~. \ivi / 

r TV/AxSn Af/A^2 ^ 

(35) 



(^^)..e = -4iVcMy^^, (33) 

n(0) 



r '^'^, , ^ = = -^ ha" - MHna' + (1 - 21n2)M2 + M^IuMH , (34) 



reg 
vac 


27r2 j 




■2a' - 


r ^ /A-\2 1 

M' N V- ' A' 


■In a' 


+ 


r ^ 




■(1 


-2 In 2) - 
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N , 


m) '^.-'"^^ 



which clearly shows that (V^V^)^ac contains quadratic and logarithmic divergences as a going to 
infinity. These divergences can respectively be removed if the subtraction constants are chosen 
to satisfy the following conditions : 

N 

M^-Y.c,k\ = 0, (36) 

i=l 

N 

M^-Y^cAt = 0. (37) 

Using the first of these conditions, the finite part of {^'4')^ac ^^"^ ^^^° ^^ expressed as 



(^^).ac = ^^Ec,f-^) logf-^) (3. 



27r2 ^ '\MJ ^\MJ 

It is now obvious that the single-subtraction Pauli-Villars scheme cannot satisfy both condi- 
tions (36) and (37) simultaneously. Although the quadratic divergence may be removed, the 
logarithmic divergence remains in {'4''4')vac ^"^^ consequently also in S{r = oo) in view of the 
relation (30). To get rid of both these divergences, we need at least two subtraction terms, 
which contains four parameters Ci, C2 and Ai, A2. The strategy for fixing these parameters is as 
follows. First by solving the two equations (36) and (37) with N = 2 for ci and C2, we obtain 



ci = (-) ^^TfAf' (39) 

C2 = -(-) T^-ry, (40) 



M\2a2_m2 
Tj AT^Af 

M\2A2_M2 

Tj Ai-A? 

which constrains the values of ci and C2, once Ai and A2 are given. For determining Ai and A2, 
we can then use two conditions (16) and (38), which amounts to adjusting the normalization 
of the pion kinetic term and the value of vacuum quark condensate. 



3 Numerical Results and Discussion 



3.1 Single- versus double-subtraction Pauli-Villars regular izat ion 

The most important parameter of the CQSM is the dynamical quark mass M, which plays the 
role of the quark-pion coupling constant thereby controlling basic soliton properties. Through- 
out the present investigation, we use the value M = 400 MeV favored from the previous 
analyses of static baryon observables. In the case of single-subtraction Pauli-Villars scheme, 
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the regulator mass Mpy is uniquely fixed to be Mpv = 570.86 MeV by using the normal- 
ization condition (12) for the pion kinetic term, and there is no other adjustable parameter 
in the model. In the case of double-subtraction Pauli-Villars scheme, we have four regular- 
ization parameters Ci,C2,Ai, and A2. From the divergence free conditions (36) and (37), ci 
and C2 are constrained as (39) and (40), while Ai and A2 are determined from (16) and (38) 
with /^ = 93 MeV and < ipip >yac — ~ (286.6 MeV) . In spite of their nonlinearity, the two 
conditions (16) and (38) are found to uniquely fix the two parameters Ai and A2 within the 
physically acceptable range of parameters. The solution that we found is 

ci = 0.445, C2 = -0.00612, Ai = 630.01 MeV, A2 = 1642.13 MeV. (41) 

As usual, all the numerical calculations are carried out by using the so-called Kahana and 
Ripka basis [17]. Following them, the plane-wave basis, introduced as a set of eigenstates of 
the free hamiltonian Hq = ol ■ V/i + /5M, is discretized by imposing an appropriate boundary 
condition for the radial wave functions at the radius D chosen to be sufficiently larger than the 
soliton size. The basis is made finite by including only those states with the momentum k as 
k < kmox- The eigenvalue problem (22) is then solved by diagonalizing the Dirac hamiltonian 
H in the above basis. We are thus able to solve the self-consistent Hartree problem and also 
to calculate any nucleon observables with full inclusion of the sea-quark degrees of freedom. 
If the theory is consistently regularized, final answers must be stable against increase of kmax 
and D (especially against the increase of kmax)- 

Now we show in Fig.l the kmax dependence of the theoretical pseudoscalar and scalar quark 
densities in the single-subtraction Pauli-Villars scheme. These curves are obtained for a fixed 
value of D as MD = 12. The corresponding kmax dependence of the quark densities in the 
double-subtraction Pauli-Villars scheme are shown in Fig. 2. Comparing the two figures, one 
immediately notices that the quark densities obtained in the single-subtraction Pauli-Villars 
scheme do not cease to increase in magnitudes as kmax increases. Undoubtedly, this must 
be a signal of logarithmic divergences contained in S{r = 00) (and generally also in P{r) 
and S{r)). On the other hand, in the case of double-subtraction Pauli-Villars scheme, the 
magnitudes of P{r) and S{r) are seen to grow much more slowly. To convince more clearly 
the above qualitative difference of the two regularization schemes, we plot in Fig. 3 the value 
of S{r = 00), i.e. the scalar quark density at the spatial infinity, as functions of kmax, and 
also as functions of log{kmax/M). Contrary to the case of single-subtraction scheme in which 
a clear signal of logarithmic divergence is observed, the value of S{r = 00) obtained in the 
double-subtraction scheme is seen to converge to some limiting value. Although the rate of 
this convergence is rather slow, it appears that this limiting value certainly coincides with the 
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prescribed value of vacuum quark condensate {4'^)yac ~ '~ (286.6 Me V)' = — 3.062 fm^ . 

Now that one has convinced the fact that the naive Pauh-Villars scheme with the single- 
subtraction term contains logarithmic divergence in the quark densities appearing in the soli- 
ton equation of motion, one may come to the following question. Why could the authors 
of ref. [12] obtain self-consistent soliton solutions despite the presence of the above-mentioned 
divergences? The answer lies in the way of obtaining a self-consistent soliton profile in the 
nonlinear model (not in the original Nambu-Jona-Lasinio model). After evaluating the pseu- 
doscalar and scalar quark densities with some (large but) finite model space (especially with 
finite kmax), a new profile function F{r) to be used in the next iterative step is obtained from 
(25). Since P{r) and S{r) appears respectively in the numerator and denominator of the argu- 
ment of arctangent, it can happen that the logarithmic divergence contained in both of P{r) 
and S{r) are offset each other. (We point out that the effect of the term f^m'^/M accompany- 
ing the scalar quark density is rather small, anyway.) In fact, Fig.4 shows the kmax dependence 
of the self-consistent profile function F{r) in both of the single-subtraction scheme and the 
double-subtraction scheme. One sees that the resultant F{r) is quite stable against the increase 
of kmax even in the single-subtraction scheme, in spite of the fact that it shows logarithmically 
divergent behavior for both of P{r) and S{r). Undoubtedly, this is the reason why the au- 
thors of [12] succeeded in obtaining self-consistent soliton profile F(r) despite the divergences 
remaining in each of P{r) and S{r). Because of this fortunate accident, self-consistent soliton 
profiles F{r) in the nonlinear model can be obtained with a good accuracy by using a modest 
value of kmax not only for the double-subtraction scheme but also for the single-subtraction 
one, and besides the resultant F{r) and not much different in these two schemes. This also 
applies to most nucleon observables which depend only on F{r) and have no direct dependence 
on S{r) and/or P{r). The previous calculation of parton distributions with use of the single- 
subtraction Pauli-Villars scheme may be justified in this sense [1-6]. To verify the validity of 
this expectation, we investigate the kmax dependence of a typical nucleon observable which 

(3) 

contains only a logarithmic divergence, i.e. the isovector axial- vector coupling constant g]^ . 

(s) 
Fig. 5 show the kmax dependence of g\ in the single- and double-subtraction Pauli-Villars reg- 

ularization schemes. One sees that this quantity certainly shows a tendency of convergence in 

both regularization schemes, though the rate of convergence in the double-subtraction scheme 

is much faster than for the scalar and pseudoscalar densities in the same regularization scheme. 

Nonetheless, one must be very careful if one is interested in nucleon observables, which have 

direct dependence on S{r) or P{r). The most important nucleon observable, which falls into 
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this category, is the nucleoli scalar charge (or the quark condensate in the nucleon) given by 

(iV|^^|A^) = / d\ [S{r) - S{r = oo)]. (42) 



The superiority of the double-subtraction scheme to the single-subtraction one must be self- 
explanatory in this case, since this quantity is convergent only in the former scheme. 



3.2 Pauli-Villars versus proper-time regular izat ion 

How to introduce ultraviolet cutoff into our effective chiral theory is a highly nontrivial problem. 
Diakonov et al. advocated the Pauli-Villars subtraction scheme as a "good" regularization 
scheme for evaluating leading-twist parton distribution functions of the nucleon within the 
chiral quark soliton model [1,2]. The reason is that it preserves several general properties of 
the parton distributions (such as positivity, factorization properties, sum rules etc.), which 
can easily be violated by a naive ultraviolet regularization. On the other hand, Schwinger's 
proper-time regularization has most frequently been used for investigating low energy nucleon 
properties within the chiral quark soliton model [7-9]. One might then wonder how these 
predictions obtained by using the proper-time regularization scheme would be altered if one 
uses the Pauli-Villars one. 

Before entering into this discussion, we think it useful to recall some basic properties of 
the proper-time regularization scheme. In this scheme, the regularized effective meson action 
takes the same form as (10) except that S^r'^^[U] is now given in the form : 



1 , . f°° dr 



S}''[U] = ^zN,J^ :^^(r)Sp(^e-^^'^-e--^oA,j^ (43) 

with 

D = i^-MU"'\ Do = i^-M. (44) 

The regularization function ^{t) is introduced so as to cut off ultraviolet divergences which 
now appear as a singularity at r = 0. For determining it, we can use a similar criterion as 
what was used in the Pauli-Villars scheme. That is, we require that the regularized theory 
reproduces the correct normalization of the pion kinetic term as well as the empirical value of 
the vacuum quark condensate. This gives two conditions : 



, , / V^(r)e^-- = f^, (45) 

^^X 7^^^^^^ = i^^)vac- (46) 
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Schwinger's original choice corresponds to taking 

V^(^) = ^ (^ - X^) ' (47) 

with A being a physical cutoff energy. However, this simplest choice cannot fulfill the two 
conditions (45) and (46) simultaneously. Then, we use here slightly more complicated form as 

V.(r) = ce(^r-^l^+{l-c)e(^r-^^, (48) 

which contains three parameters c, Ai and A2 [18]. Although the above two conditions are 
not enough to uniquely fix the above three parameters, we find that solution sets (c, Ai,A2) 
lie only in a small range of parameter space and that this slight difference of regularization 
parameters hardly affects the soliton properties. We use the following set of parameters in the 
numerical investigation below : 

c = 0.720, Ai = 412.79 Me V, A2 = 1330.60 MeV. (49) 

Within the framework of the chiral quark soliton model, which assumes slow collective rotation 
of a hedgehog soliton as 

U"'^^{x,t) = A{t)U^"{x)A\t), A{t) c SU(2), (50) 

the nucleon matrix element of any quark bilinear operator i/jOiIj is given as a perturbative 
series in the collective angular velocity operator Q defined by 

n = iA\t)jA{t). (51) 

It is shown below that a noteworthy difference between the proper-time regularization and the 
Pauli-Villars one appears at the zeroth order term in Q. We recall that, in both schemes, the 
0{fi^) contribution to this matrix element is given as 

{Of = JVA ^£mJ^] (O) J ^SmJ^]' (52) 



with 



' vp 

,(J) 



{0f = {0C + {0)Z^ (53) 



where "^ m m t^] ^^ ^ wave function describing the collective rotational motion. In eq.(53), 

{O)Tai = Nc{0\d\0), with d = A^OA, (54) 
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represents the contribution of the discrete bound state level called the valence quark one. 
Within the Pauli-Villars scheme, the contribution of the Dirac continuum can be given in 
either of the following two forms : 

(^)i;p ~ ^c X] (^1^1^) ~ Pauli-Villars subtraction, 

?i<0 

= — Nc ^^ {n\0\n) — Pauli-Villars subtraction. (55) 

n>0 

Note that the first form is given as a sum over the occupied single-quark levels, while the second 
form given as a sum over the nonoccupied levels. The equivalence of the two expressions follows 
from the identity 

= Spd = J2{^\O\n) + J2{^\O\n), (56) 

n<0 n>0 

which holds for most operators including the isovector magnetic moment operator investigated 
below, if it is combined with the fact that a similar identity holds also for the corresponding 
Pauli-Villars subtraction terms. The situation is a little different for the proper-time regular- 
ization scheme. The regularized Dirac sea contribution in this scheme is given in the following 
form [8] : 

{OC = -Y E sign(E„)^(i?„)(n|0|n), (57) 

n=all 



1 f^ dr 



2 



g{En) = ^ ^\En\e-^^-. (5J 

a/tt Jo Jt 



with 



To compare this with the corresponding expression in the Pauli-Villars scheme, it is convenient 
to rewrite it as 

(OC = \{NcY.9i^n){n\6\n)-N,Y.9i^n){n\b\n)\. (59) 

^ n<0 n>0 

One sees that here the answer is given as an average of the two expressions, i.e. the one given 
as a sum over the occupied levels and the others given as a sum over the nonoccupied levels. 
(This feature is a consequence of the starting covariant expression for an operator expectation 
value in the proper-time scheme.) However, contrary to the previous case in which ultraviolet 
regularization is introduced in the form of the Pauli-Villars subtraction, now there is no reason 
to believe that the above two terms give the same answer. In fact, the introduction of the 
energy dependent cutoff factor giEn) generally breaks the equivalence of the two expressions 
because of the spectral asymmetry of the positive- and negative-energy levels induced by the 
background pion field of hedgehog form. 
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Now we start a comparative analysis of the two regularization schemes on the basis of the 
numerical results. For reference, we also solve the soliton equation of motion in the chiral limit. 
By assuming no (or at least weak) m.„ dependence of < ipi/j >yac appearing in (16) and (38), 
this calculation can be done by setting m^ = in (18) and (25) without changing the sets of 
regularization parameters given in (41) and (49). Since the way of cutting off the ultraviolet 
component is totally different for the two regularization schemes, it naturally affects solutions 
of the soliton equation of motion. Although the detailed contents of the soliton energy are 
highly model dependent concepts and are not direct observables, they are anyhow very sensitive 
to this difference of the self-consistent solutions. Table 1 shows this comparison. Comparing 
the answers of the two regularization schemes, one finds that the Pauli-Villars scheme leads 
to more strongly deformed soliton, which means a deeper binding of the discrete valence level 
and larger vacuum polarization energy. One sees that the total soliton energy is lower for 
the Pauli-Villars scheme than for the proper-time scheme. One also observes that the soliton 
energy is very sensitive to the pion mass. When one goes from the finite pion mass case to the 
chiral limit, one obtains much lower soliton energy. 



Table 1: The static soliton energy in the proper-time regularization schme and the (double- 
subtraction) Pauli-Villars one. E^ai, El^^ respectively stand for the valence quark contribution 
and the Dirac sea one to the fermionic energy, while E^ represents the mesonic part of the 
energy. The sum of these three parts gives the total static energy El^^^uc- 





E^ai [MeV] 


EIJ [MeV] 


E„, [MeV] 


KtLc [MeV] 


proper-time (?72^ = 138 MeV) 
Pauli-Villars (m^ = 138 MeV) 


633.0 

447.6 


617.6 
569.2 




37.2 
51.3 


1287.9 
1068.1 


proper-time (m^r = MeV) 
Pauli-Villars (m^ = MeV) 


555.6 
351.5 


688.6 
655.4 








1244.2 
1006.9 



Probably, the most important observable which has strong sensitivity to the above difference 
of the self-consistent solutions is the flavor-singlet axial charge or the quark spin content of the 
nucleon (S3). The theoretical predictions for this quantities in the two regularization schemes 
are shown in Table 2. In evaluating this quantity, we did not introduce any regularization, 
because it is related to the imaginary part of the (Euclidean) effective action and is convergent 
itself. This means that the difference between the two schemes purely comes from that of 
the self-consistent solutions. One sees that the Pauli-Villars scheme leads to smaller quark 
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Table 2: The quark spin content of the nucleon < S3 > in the proper-time regularization 
scheme and the Pauh-Villars one. 







< S3 >^al 


< S3 >v.p. 


<S3> 


proper-time {771^^ = 
Pauh-Villars {rrij, = 


: 138 MeV) 
= 138 MeV) 


0.484 
0.391 


0.005 
0.008 


0.489 
0.399 


proper-time (m^ -- 
Pauli-Villars (m^ 


= MeV) 
= MeV) 


0.374 
0.286 


0.007 
0.011 


0.380 
0.298 



spin content. The reason can easily be understood. Within the framework of the chiral quark 
soliton model, the rest of the nucleon spin is carried by the orbital angular momentum of quark 
fields and this latter portion increases as the deformation of the soliton becomes larger [8]. A 
similar tendency is also observed when one goes from the finite pion mass case to the chiral 
limit. 

Table 3: The 0{Q^) contributions to the isovector magnetic moment of the nucleon in the 
proper-time regularization scheme and the Pauli-Villars one. The second column represents 
for the valence quark contribution. The third and fourth columns stand for the answers for 
the vacuum polarization contributions respectively obtained with the occupied and nonoccu- 
pied formulas, while the fifth column gives the average of the two answers. The total 0{Q^) 
contributions are shown in the sixth column. 







f^^lm 


occupid 


nonoccupied 


average 


/i(3)(l]0) 


proper-time {rrij^ = 
Pauli-Villars {rrij^ = 


■ 138 MeV) 
= 138 MeV) 


1.611 
1.762 


1.312 
0.996 


0.210 
0.996 


0.761 
0.996 


2.372 
2.759 


proper-time (m,r '■ 
Pauli-Villars {rrij^ 


= MeV) 
= MeV) 


1.623 
1.810 


1.908 
1.738 


0.588 
1.738 


1.248 
1.738 


2.875 
3.547 



There are different kinds of nucleon observables, which contain (potential) logarithmic 

divergence and thus depend directly on how they are regularized. Most typical are the 0{Q^) 

(3) 
contribution to the isovector axial- vector coupling constant g)^ and the isovector magnetic 

moment /i*^^-* of the nucleon. Let us first show the results for the isovector magnetic moment, 
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since it turns out to have stronger dependence on the choice of the regularization scheme. Table 
3 shows the 0{Q^) contribution to the isovector magnetic moment. For each regularization 
scheme, the second column represents the answer obtained with the occupied expression, while 
the third column does the answer obtained with the nonoccupied one. In the case of Pauli- 
Villars scheme, the equivalence of the two expressions is nicely confirmed by the explicit 
numerical calculation. In the case of proper-time scheme, however, we encounter quite a 
dissimilar situation. First, the answer obtained with the occupied expression is about 30 % 
larger than the corresponding answer of the Pauli-Villars scheme, while the answer obtained 
with the nonoccupied expression is about 80 % smaller than the answer obtained with the 
occupied one. Since the final answer of the proper-time scheme is given as an average of the 
occupied and nonoccupied expressions, the consequence is that the prediction of the proper- 
time scheme for the 0{rf) contribution to jj,'^^^ is about 14% smaller than the corresponding 
prediction of the Pauli-Villars scheme. (See the fourth column of the Table 3.) Note that the 
difference between the two regularization schemes becomes much more drastic when one goes 
to the chiral limit. This is due to the fact that the 0{Q^) vacuum polarization contribution 
to the isovector magnetic moment is extremely sensitive to the pion mass effect such that it is 
much larger in the chiral limit. 



Table 4: The final predictions for the isovector magnetic moment of the nucleon, given as sums 
of the 0{fi^) and 0{fl^) contributions. 







/i(3)(f]0) 


/i(3)(f]i) 


/i(3)(fi0 + l]l) 


proper-time (m^ = 
Pauli-Villars (m^ = 


: 138 MeV) 
= 138 MeV) 


2.372 

2.759 


1.072 
1.211 


3.445 
3.970 


proper-time (m,r '■ 
Pauli-Villars (m^ 


= MeV) 
= MeV) 


2.875 
3.547 


1.032 

1.182 


3.907 

4.729 



Before comparing our theoretical predictions with the observed isovector magnetic moment 
of the nucleon, we must take account of the 0{Q^) contribution, too, since it is known to give 
sizable correction to the leading-order result [19,20]. Although we do not go into the detail 
here, it turns out that this 0{Q^) piece is not so sensitive to the difference of the regularization 
scheme as the 0{Q^) piece is. The reason is that this 0{Q^) term is given as a double sum 
over the occupied levels and the nonoccupied ones and the formula has some symmetry under 
the exchange of these two types of single-quark orbitals [21]. The final predictions for the 
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Table 5: The final predictions for the isovector axial-coupling constant of the nucleon, given 
as sums of the 0{Q^) and 0{Q^) contributions. 







gfm 


gfm 


(3~ 


\n^ + n^) 


proper-time (m^ = 
Pauli-Villars (m^ = 


■ 138 MeV) 
= 138 MeV) 


0.848 
0.976 


0.412 
0.408 




1.260 
1.384 


proper-time {171^^ - 
Pauli-Villars {rriTr 


= MeV) 
= MeV) 


0.921 
1.054 


0.348 
0.344 




1.269 
1.398 



nucleon isovector magnetic moment obtained as a sum of the 0{Q^) and 0{Q^) contributions 
are shown in Table 4. After all, the prediction of the Pauli-Villars scheme is about 15 % larger 
than that of the proper-time scheme and a little closer to the observed moment. The effect 
is much more drastic in the chiral limit. The prediction of the Pauli-Villars scheme is about 
20 % larger than that of the proper-time scheme and nearly reproduces the observed isovector 
magnetic moment of the nucleon, i.e. ^^J^Jp ~ 4.71. 

Finally, we show in Table 5 the predictions for the isovector axial-charge of the nucleon 
obtained as a sum of the 0{Q^) and 0{Q^) contributions. Also for this quantity, there are some 
detailed differences between the predictions of the two regularization schemes. Nonetheless, 
the final answers for g]^ turn out to be not so sensitive to the difference of the regularization 
schemes as compared with the case of the isovector magnetic moment. Besides, one also notices 
that the finite pion mass effect hardly influences the flnal prediction for this particular quantity. 



4 Conclusion 

In summary, the single-subtraction Pauli-Villars regularization scheme, which is often used in 
evaluating nucleon structure functions within the framework of the CQSM, cannot be regarded 
as a fully consistent regularization scheme in that it still contains ultraviolet divergences in the 
scalar and psuedoscalar quark densities appearing in the soliton equation of motion. However, 
these divergences can easily be removed by increasing the number of subtraction term from one 
to two. After this straightforward generalization, the effective theory is totally divergence free. 
Especially, both the vacuum quark condensate and the isoscalar piece of the nucleon scalar 
charge becomes flnite now. Nonetheless, we flnd that, owing to the accidental cancellation 
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explained in the text, one can obtain a finite soliton profile F{r) even in the single-subtraction 
scheme, and besides the resultant soliton solution is not extreme different from the correspond- 
ing one obtained in the double-subtraction scheme. Furthermore, it turns out that, for most 
nucleon observables, which contain only the logarithmic divergence, the predictions of the two 
regularization schemes are not much different. The previous calculations of quark distribu- 
tion functions with use of the single-subtraction Pauli-Villars regularization scheme would be 
justified in this sense. 

We have also carried out a comparative analysis of typical nucleon observables based on 
the Pauli-Villars regularization scheme and the proper-time one. A nice property of the Pauli- 
Villars regularization scheme, which is not possessed by the proper-time one, is that it preserves 
a nontrivial symmetry of the original theory, i.e. the equivalence of the occupied and nonoccu- 
pied expressions for 0{Q^) contributions to nucleon observables. The improvement obtained 
for the isovector magnetic moment of the nucleon was shown to be related to this favorable 
property of the Pauli-Villars regularization scheme. How to introduce ultraviolet cutoff into 
an effective low energy model should in principle be predictable from the underlying QCD 
dynamics. For lack of precise information about it, however, phenomenology must provides us 
with an important criterion for selecting regularization schemes. The regularization scheme 
based on the Pauli-Villars subtraction appears to be a good candidate also in this respect. 
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Figure caption 



Fig. 1. The kmax dependence of the scalar quark density S{r) and the pseu- 

doscalar density P{r) in the single-subtraction Pauli-Villars scheme. 

Fig. 2. The kmax dependence of the scalar quark density S{r) and the pseu- 

doscalar density P{r) in the double-subtraction Pauli-Villars scheme. 

Fig. 3. The scalar quark densities at the spatial infinity S{r = oo) as functions 

of kmax/M and as functions of \og{kmax/M) in the single- and double- 
subtraction Pauli-Villars schemes. 

Fig. 4. The kmax dependence of the self-consistent soliton profiles F{r) in the 

single- and double-subtraction Pauli-Villars schemes. The curves with 
different kmax are almost indistinguishable. 



(3) 

Fig. 5. The kmax dependence of the nucleon isovector axial-charges g]^ in the 

single- and double-subtraction Pauli-Villars schemes. 
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